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Abstract 

We consider a quantum linear oscillator coupled to a bath in equilibrium at an arbitrary temperature and then exposed 
to an external field arbitrary in form and strength. We then derive the reduced density operator in closed form of the 
coupled oscillator in a non-equilibrium state at an arbitrary time. 
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1. Introduction 

One of the most successful approaches to non-equilibrium statistical mechanics is the linear response theory 
H] 121111. This allows us to predict the average response of a physical quantity of the system to external perturbations 
with weak strength. At the heart of hnear response theory we have the fluctuation-dissipation theorem [4. J |, which 
offers a clear-cut relationship between irreversible processes in a non-equilibrium state and thermal fluctuations in the 
(initial) equilibrium state. However, this approach is, in general, restricted to non-equilibrium states near equilibrium 
in validity. 

As a prototype of quantum dissipative systems, the scheme of quantum Brownian motion has been studied deeply 
and widely over a long period f5\. At its heart we have a quantum harmonic oscillator linearly coupled to an 
independent-oscillator model of a heat bath [quantum Brownian oscillator] in equilibrium at a (low) temperature. 
Due to its mathematical simplicity, this system allows the linear response theory to yield an exact expression for an 
average response of the system operator q (position) [and also that of p (momentum)] to external forces F arbitrary in 
form and strength as well as those for the equilibrium fluctuations and {p^)i3, respectively ||5]|6|. Based on this 
well-known result, quantum Brownian oscillator has recently attracted considerable interest in investigating thermo- 
dynamic behaviors of small-scaled quantum systems coupled to quantum environments in the low-temperature regime 
Q [8. 9^^. 1& JT] [12] [T3l 14| ("quantum thermodynamics" lfT5l[T6llT7J ). Here the finite coupling strength between sys- 
tem and environment yields some quantum subtleties and so cannot be neglected whereas ordinary quantum statistical 
mechanics is intrinsically based on a vanishingly small coupling between them. 

On the other hand, response functions in a far-from-equilibrium state such as {q"{t)) and {p"(t)} with n > 2 in this 
system cannot be obtained directly from the linear response theory. To explicitly have such non-equilibrium quantities, 
we need to exactly treat the higher-order terms in the external field and accordingly go beyond the scheme of linear 
response theory [18, .19J . The primary goal of this paper is to derive a reduced density operator in closed form of 
the coupled oscillator in a non-equilibrium state at an arbitrary time f, which can, obviously, provide all higher-order 
fluctuations of the non-equilibrium state. For numerical analysis we will also consider a variety of external fields (d.c. 
and a.c.) leading to explicit evaluation of the non-equilibrium fluctuations. In doing so, we will employ not only the 
qF(t) interaction Hamiltonian ("scalar-potential gauge") but also the pA{t) interaction Hamiltonian ("vector-potential 
gauge"). The equivalence of the two interactions is based on the gauge transformation between their wavefunctions 
satisfying the corresponding (time-dependent) Schrodinger equations, respectively (for a detailed discussion of qF 
versus pA gauge problem, see Ref. EOl ). And we will adopt the Drude model with a finite frequency cut-off' for the 
spectral density of bath modes, which is a prototype for physically realistic damping ||5l. 

The general layout of this paper is the following. In Sect. |2] we briefly review the general results of quantum 
Brownian oscillator needed for our later discussions. In Sect. [3]we explicitly derive an exact expression for the reduced 
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density operator of the coupled oscillator in a non-equilibrium state. In Sect. |4]we perform numerical analysis for 
various non-equilibrium quantities induced from the reduced density operator Finally we give the concluding remarks 
of this paper in Sect. |5] 



2. General treatment of quantum Brownian oscillator 

The quantum Brownian oscillator under consideration is described by the model Hamiltonian ||5]|6l 

H - Hs + Hh-sb , 



where 
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The Hamiltonian Hssb can split into the bath and the coupling terms such as 
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respectively. The interaction Hamiltonian Hsb contains a term of quadrature completion modifying the frequency 
of the system oscillator Hs- This quadrature term is needed not only for later exact calculations; in fact, the model 
Hamiltonian without the quadrature completion has a significant defect, as was discussed in IISTI . The total system is 
assumed to be within the canonical thermal equilibrium state pp - e^l^^ jZp where /? = lliks T), and Zg is the partition 
function. 

From the Heisenberg equations of motion for q and p we can derive the quantum Langevin equation |5, 6| 



M§(f) + Mj^dTy(t-T)^(T) + Mojlq(t) = 1(f), 
where we used /3(f) - Mqit), and the damping kernel and the noise operator are, respectively, given as 
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Here the expectation value of the noise operator vanishes, Tr{f(f)p^) - or, equivalently, (|(f))p,y - with re- 
spect to the initial bath state, prepared as a shifted canonical equilibrium distribution, pb' = g-Z'f * * f^^' ) jjj which a 
normalization constant zj^ ^ is the properly defined partition function [SJ. And the noise correlation is given as 181 

1 h ^ ' 

S^dt - O 2 + ^(^')^(^))p„ = 2 ^ ''"'''^^■^^ " coth(^) . (8) 

With ^ — > 0, the correlation S^^t - f) reduces to its classical counterpart, My{t - t')//3 [22J. We also introduce a 
response function ||6l 

WO = i <m§(o)])/,0(f), (9) 

where 0(f) represents a step function. Then we can have other response functions as well such as Xpgi^) - ~Xgp(f) - 
Mxqqit), and;t'„,(f) = -M^Xqq(f). 



For a later purpose it is instructive to discuss the time-reversal dynamics of q{t) in terms of r(t) :- q(-t) and its 
momentum s{t) :- -p(-t). We can then derive the corresponding quantum Langevin equation fS) 



M ht) + M dTjit - r) Ht) + McjI f(t) = |-(0 ■ (10) 

While this is the same in form as equation (j6|, the two equations differ in the noise in such a way that |-(r) is identical 
to §it), however, with the replacement of Pj(0) — » -pj(0) in (j7|. And from equation ^ and stationarity of the 
equilibrium correlation function between operators F and G such as {F(t) G(0))/j - {F(Q) G{-t)}j3 f^, we can easily 
obtain ;^f,.,.(f) = -Xqq(t) [note that r(0) = qiO)]. Likewise, it also appears that Xrp(t) -XqpiO- 

Now we intend to derive explicit expressions for q(t) and r(f), respectively. To do so, we first apply the Laplace 



transforms to equations (j6j) and (10 1, respectively. Let the Laplace transform q(s) :- £.{q(t)}is), and so we have 
X,{q{t)}(s) = sq{s) - q{0) and X.{q{t)}{s) - s^qis) - sq(0) - q{Q) [23J. We can then obtain the Fourier-Laplace 
transform of q{t), which reads as 

q(oj) ^ , 2 q - h ^ ■ ~. . ^ P - 11 2 ■ I . y ii<^Xi - ^'l (11) 

i 

where q{u)) - q(s) with s - -io) + 0^, and the Fourier-Laplace transform of r(f), 

r(w) = . . , , q + , . ' . p - h . ■ I . ^ y -T^ ('wJc/ + ^\ (12) 
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where £((y) - l{s). The operators q,p,Xj, and pj represent the initial values q{0), p{Q), Xj{Q), and Pj{Q), respectively. 
And the Fourier-Laplace transform of y(f) is 

= (13) 

] ' 1 J 

We introduce the susceptibility, defined as the Fourier-Laplace transform of Xqq{t) in (j9]l, such as |l5]|6l 

Xqq{(^) := fl'f;r«?(Oe'"' = {mu),q])p, (14) 
which easily reduces to -l/{M(w^ + icofiai) - Wy)) with the aid of ( 1 1 1. And it then appears that 

Xq.yi<^) = ii{[q((^),Xj]}fi = ,„.j2_^2) ^w(^)' (15) 
Xrriui) = ^ <[£(w),§])/j = -Xqqi^)^ (16) 

Xqpi^^) = iojMxqq(co), and Xqpj((^) = iajmjXqxj{(^)- Likewise, we also have Xrpit^) = Xqp(<^)^ Xrxj((^) = -jr?x/w), and 
Xrpj{(^) - Xqpji^)- Consequently, equations (111 and ( 12 1 can be rewritten in terms of the susceptibilities, respectively, 
as the compact expressions 

J?(W) = -Xqpi'^)q+Xqqi'^)p-^\Xqp,i(^)^i-Xqx,{(^)Pj\ (17a) 
Ki^) = -^r;Xw)§-Hi-„.(w)^-^{^rp,(w)x— ^„,(t^)^;)- (17b) 

Applying the inverse Fourier transforms to \ lla\ and ( 17b i immediately allow us to have 

q(t) = -Xqp(t)q+Xqqit)p-^{XqPj(t)^j-XqXi(t)pj] (18a) 

j 

rit) = -Xrp(f)q+Xrr(f)p-Yp:'-P^'^^j-Xrx,{t)pj], (18b) 

j 
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where 

Xmit) = Tn n dcoxmioj) e'^^' . (19) 

The two equations will be used in Sect. |3]for derivation of the reduced density operator of the coupled oscillator in 
closed form. 

For a later purpose it is useful to introduce well-known expressions for the equilibrium fluctuations in terms of the 
susceptibility ^^^(w) such as ll24ll 

<§'>-8 = f r (^) lm\x,,{cn + i 0+)) (20) 

and 

<P')/S = ^ r '^o* (^) lm{x,,(oj + i 0+)) , (21) 

respectively, which can be derived from the fluctuation-dissipation theorem HIS]. 



3. Reduced density operator of the coupled oscillator in a non-equilibrium state 

Now we study the influence of an external field on a linear oscillator coupled to a bath. To do so, we first consider 
the equation of motion for the density operator of the total system {i.e., oscillator plus bath), which reads |l5]|6| 

p(f) = e-"^op(0) - iJ^dTe-'^'-^^^o£i{T)p(T). (22) 



Jo 

In the scalar-potential gauge for the field-coupling, the total Hamiltonian reads as 'Hs(t) - ft - qF(t), and the corre- 

-I 

£op(t) = ^[H,p(t)]; 1\'\t)p(t) = -l[q,p(T)]F(T). (23) 



sponding Liouville operator X = Xo + X*,'^ satisfies 



In the vector-potential gauge, on the other hand, the total Hamiltonian is given as 

= + 1 + ^'-'i' ' ^^^^ 

which is identical H + pdt) p/M + p^{t)/(2M), and accordingly 

to Pit) = i[H,p(t)] ; Xf*(T)p(T) = ^[/3,p(t)] ^ . (25) 

Here the vector potential pdt) = J^^ F(T)dT is an impulse induced by field F. For an uncoupled oscillator the 
equivalence of the two interactions is based on the gauge transformation ij/siq,!) - e^^ A*')^ f2^, where the 
wavefunctions t//s(q, t) and i/'j(^, t) satisfy the time-dependent Schrodinger equations in the scalar-potential and the 
corresponding vector-potential gauges, respectively. 

We substitute (23 1 into (22 1, with p(0) = pp, and make iterations forp(T) in the integral. Then we can arrive at the 
expression 

p,(f) - PP+ i /o dTF{T) [q,pp\ ei^'-'^" + [if /J drFir) £ dr' F(t') e-rJ'-''>" X 

[§,e-s(^-^')H[§^P^jei(r-T')Hjei(/-r)H + (1)^ JJ^fTF(T) £ dx' Fir') £' dr" F(T")e-i^'-^^" X 

[§,p^]e^(-'-")H]ei(T-T')Hjei('-T)H + . . . . (26) 

With the aid of [pg, H] - 0, this equation easily reduces to the expression in terms of r(f) = e^i'^qei'^ such as 

pAt) ^ PP+ { /o dT F{t) [Ht -T),pp\+{{f X' dr Fir) £ dr' Fir') [r(f - t), [r(t -T'),pp]] + 

[jrf £ dTF(T) £ dT' F(t') £' dT" F(t") [r(t - r), [r(f - r'), [r(f - r"),P/,]]] + ■ ■ ■ , (27) 
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where 

[r{t),Pfi] = -Arr;,(f)[§,P/;] + Xrr(t)[p,Pl}] ' J^XrpjiO [xj, pfl] + ^A-rx/O [^j^P/s] , (28) 

J j 

obtained directly from equation ( |18b| l. For the vector-potential gauge, on the other hand, we plug ( [25] l into ( [22] l, 
and after making some calculations similar to those for (|26]) we can finally obtain the expression in terms of s(t) - 
-g-ji'^ p g-k'ti such as 



From (28 i and s(t) = Mr(f) we can easily have 

[s(t),P/}] = -A'sp(f)[§,P/s] + A'.5r(f)[p,P/5] - X-¥.5p/0[-^j,P/j] + ^A's.v/Of^j-P/?] ■ (30) 



Let us now consider the reduced density operators for the coupled oscillator, ??(?) := Tti, p{t) from (27 1 and (29 1, 
respectively. Here, Tr^ denotes the partial trace for the bath alone. The initial state '^(0) of the coupled oscillator, 
being the reduced operator of the canonical equilibrium state pp, is known as ||5l l25l 



First, from (|28l) we can obtain 



" {q\Txk [mrpp] W) = S,,.{t) {q\m)\q') , (32) 

where Sqq'{t) :- -ihxrrit) {dq + dqi^ - Xrp(t) il - and similarly from (30i we can also have 

(qlTTb [m,p4 W) = %?'(0 {q\m)\q'} , (33) 

where %^'(f) := -ihxAt) {dq + dq') -Xspit) (q - q')- Here we used 

k 

rooY['^Pk{Pk\[PhPp\\Pk) = 0. (34b) 

k 

From ( |18b[ ), ( |34a| i and \3Ah\ it also appears that 

n n ^-"^ |[^' [^'^pp\\ = (35) 

unless both (3 e {q,p} and (3' € {q,p}. Therefore we can arrive at the expressions, 

<^|Tr4r(f),[r(T),p;5]]|^') = Sqq,(t)Sqq,(T){q\m)\q') (36) 

<9|Tr4i(f),[.?(T),p;5]]|^'> = 'Vqq,(t)^qq,(T){qm)\q'). (37) 



Here we also used 



{q\Trh [p, [p,pp]] \q') = (f )' (dq + dq-f {qm)\q'} 
{qlTri, [q, [p,p^]] \q') ^ (q - q') [dq + dq,) {q\m)W) 

{q\TXh [q, [q,pp^ W) ^ (q - q'f (qimW) , (38) 
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and [q, [p,Pfi]] = [p, [q,p/3]]- Along the same line, after making lengthy calculations, we can also obtain 

<<?|Tr,[r(f),[r(T),[r(T'),P/,]]]|^'> = Sg,it)S,g,(T)S,g,(T'){qm)\q'} (39) 
{q\Tvb [s(tl [s(t), [s(t'),PiI]]] \q'} = %,-(f) %,-(t) %,-(t') <^|iR(0)|^') . (40) 



With the help of equations (27 1, ( |32) i, (36 1, and ( [39| we can finally find the matrix elements of the reduced density 
operator in the scalar-potential gauge such as 

{q\kMq'} = T e^'''''Hq\'k(mq'} , (41) 

where the operator Js(t) dr F(t) Sqq'(t-T) represents a time-evolution action, and T is the time ordering operator. 
For the vector-potential gauge, along the same line, from ( [29| , ( [33] l, ( [37] ), and ( |40| we can arrive at the matrix elements 

(qlkm'} = T e^'' '^''> {qm)\q'} , (42) 

where the operator /,,(f) :- (1/M) dr Pc{T)'Vqq>{t - r). 

Now we simplify the above expressions for {q\'Rs{t)\q') and {qyR.i,{t)\q'). The operator 75(f) immediately reduces 
to ih{q{t))s {dq + dqi) + {p{t))s (q - q') from the well-known exact expression 

0(f)). - 0(0)) = dTF(T)xog(t - T) , (43) 

where O e {q, p}, obtained directly from the linear response theory ||6| (note here that (§(0)) - = and {p(Q)) - 



{p)l3 = ISHUl, and {p{t))s = M<§(f)),). Due to the fact that [q-q',dq+ dq'] = 0, equation (41 



then becomes 



{qrnrn') = e^W)). (,-,') £ tMML + {q\m)\q') ■ (44) 



n=0 



Let y :- q + q' so that dq + dq' = 2dy. Then we can easily obtain 

[dq . dq')" (qmOW) = r ^ exp(-^ - , (45) 

which can subsequently be expressed in terms of the Hermite polynomial as 



e-l^^i-'i"^\--- H„{z) (46) 



with z — yl -\/8{q^)jj. Here we used the identity H„(z) — (-1)" (djdzf e " ||26i . Then, with the aid of the identity 
e^~'^'^ - lijJLo{//n(z)/«!} Il26l . equation (44i finally reduces to the exact expression 



The normalization Tr:Rs(f) = 1 can easily be shown with the aid of Il26ll 

^"^dqe-^f"r^^'"'^ = yf^e''''". (48) 



For the vector-potential gauge, along the same line, after making some calculations leading to (44 1 with 7y(f) in place 
of 7.5(f), we can obtain 

00 

(qimW) = e^^P'-'^>" ^ (5^^ + dq')" {qm)\q') , (49) 

«=() 
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where we identified {q(t))y and {p(t))v, respectively, as dr Pc(T)Xqq(t-T) and J^'^ dr PcijJXpqit-T) [note ihsAXsrif) - 
—Mxqqif) wdxsp(t) - -Mxpqif)}, Verified directly from the linear response theory with (25 i in place of (23 1. In fact, 
with the aid of integration by parts it can easily be shown that 



Consequently we can immediately obtain 

<^|^,(f)|<?'> = {qmt)\q' 



e 



-T;Pc{t){q-q') 



{q\'RMq')- 



(50) 



(51) 



Obviously, we have Tr!Rv(0 = 1 from (47 1. 

It is instructive now to consider (q^it)) and (p^it)) in both scalar-potential and the vector-potential gauges. Using 
( 47 1 we can easily obtain 



<§2(f)), = {f)^ + {q(t))l 
{p\t)), = {p\ + {pit))] 



Similarly, we can also have 



{q\t)),. = Tr {§2^,(0} 

{p\t))y = {p\ + {p(t))l 



{q\t)). 



As a result, the instantaneous internal energy in the scalar-potential gauge 

is not necessarily identical to its vector-potential gauge counterpart, namely, 

<H,(f)), - <H,(f)).. = ^ [{pit))] - {pit))] 



+ 0. 



(52a) 
(52b) 



(53a) 
(53b) 



(54) 



(55) 



At first glance, it looks like a paradox. However, we have a rather simple justification for this EtII : In the scalar- 
potential gauge problem, the experiment is performed in such a way that we turn on an external field at f = and 
then turn ofi" at f = f/. Afterwards we measure the fluctuation {p^itj))s. In the vector-potential-gauge setting, on the 
other hand, we need to turn off the vector potential pAtf) - J,'^ Fir) dr rather than the external field. Consequently 
the fluctuation {p^itf))v differs from its scalar-potential counterpart in such a way that 



{p\tf)),. = {p\ + {{pitf))s-pcitf)}\ 



(56) 



which is actually accordance with the result in ( |53b| l with ( [50) l. It is, however, in general physically unrealistic to carry 
out an experiment in which the vector potential is turned off. 

Comments deserve here. First, it is interesting to note a time-independent behavior of the purity measure 



TiK^it) = TT'R]it) = 



(57) 



obtained directly from (47 1 and (51 1 with (48 i, respectively. 

Secondly, as was shortly pointed out in Sect. [11 whereas the equilibrium quantities and (/3^)/s in equations 



( [So] ) and ( [2T| can exactly be obtained from the scheme of linear response theory, this is not the case for their non- 
equilibrium counterparts {q^it)) and {p^it)); in fact, by using (OTJi we can arrive at the expression 



{qHm, = <§2)^ + j;' dTX^'Xt - T) FiT) + £ dT FiT) dr' FiT')x^^\t - r, f - r') , 



(58) 
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wheTe x^^\t) - Xq^qi^) - ([§^(0^§])/3 = 0, and the 2nd-order response function ;^'^^^(f,T) = r(f)], ?^(t)])^ can be 
obtained from the cycHc invariance of the trace, which subsequently reduces to 2xqq{f)Xqq(j) with the aid of ( 18b| l. 
The relation J^^' dr dr' - J^^ dr' J^' dr 1281 then allows equation ( 58 1 to become 



which is also identical to 



(q\t)),-{q\ = 2 J^dT'F(T')Xqq(t-T') £, dz F(t) Xqq(t - t) ., 

2 /; dr' F{r')xqq{t - r') ^ dr F{r)xqq{t - r) 



(59) 



(60) 



directly resulting from (58 i with exchange of the two integral variables t and t' . From ( |43] l, (59 1 and (60 1 we can 
immediately recover the exact result in (52ai. Similarly we can do the same job for {p^{t))s and then for their vector- 
potential gauge counterparts, respectively. 

Subsequently, we can also obtain, with the aid of ( [48] ), the higher-order fluctuations such as 

(61a) 
(61b) 



and 



{p\t)},, = 3{p^}l + 6{p^)p{p{t))l + {pit))], 
etc. Their vector-potential counterparts immediately appear with the replacement of {p(t))s — > {p{t))y. 



(62a) 
(62b) 



4. Numerical Analysis within the Drude damping model 

We carry out the numerical analysis in the scheme of the well-known Drude model (with a cut-off frequency ojd 
and a damping parameter y,,), which is a prototype for physically realistic damping. It is then known that ifTTI 

<§T - (63a) 

/=i 

<pX - -^i^vM^i-^^^C-S')} ^^^^^ 

/=1 

in terms of the digamma function tf/iy) - d \nY{y)ldy ll26ll . where lox - Q., - Z\, - Zi, and the coefficients 

i(i) _ ;i+z2 . ,(2) _ n + z2 . i(3) _ n + zi /--Hx 

Here we have employed, in place of (ljq, wj, Jo), the parameters (wo, Q, y) through the relations Q 

0.^= ^ Q + r; r. - T ^^^^^ ' ^^5) 

and then zi = y/2 + iwi and zi - 7/2 - iwi with Wi = -y/(wo)^ - (y/2)^. For the underdamped case (wq > y/2) we 
have Z2 - Z\ whereas z\, Zi > for the overdamped case (wo < y/2). The susceptibility in the Drude damping model 
is also well-known as 191 [TOl 

-Id), ^ 1 OJ + i(n + Zl + Zl) ,,,, 

xfJi^o) = - — - — — — — -— - — -— . (66) 



With the aid of ( [T9| ) we can easily obtain the response function 

X^qq(t) = - 77 • (67) 



'"'"'^ " ~M (Q-zi) (Z1-Z2) (Z2-n) 
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This is real-valued and holds true for both underdamped and overdamped cases. 

We first apply a static external field Fi (f) = A] 0(f) (d.c. field) and then an oscillatory field F2(0 = A2 sinw/f(a.c. 



field). By substituting equation (|67|i to (|43j) we can easily obtain both {q{t)}f^ and {p(t)}f^ within the Drude damping 



model in closed form for the two different forms of external force, respectively. Similarly we can also have their 
vector-potential counterparts in closed form. Figures l]|4 demonstrate temporal behaviors of (§"(0)*/' with n - 1,2 
for different damping and control parameters. Further, it is instructive to study a temporal behavior of a distance 
between the initial equilibrium state !R(0) and the non-equilibrium state fl{t). To do so, we adopt a well-defined 
measure D^(t) = Tr({:R(f) -!R(0))^), introduced in ll29ll . which is, independent of the dimension of the Liouville space, 
between and 2. With the aid of (|48l) we can then have 



(68) 



and Z)J(f) = £'?(OI(;)(o>,^(/5(0>w- figurejsjthis measure within the Drude damping model is demonstrated for different 
external fields and temperatures. 



5. Concluding remarks 

In summary, we have discussed the field-induced dynamics in the scheme of quantum Brownian oscillator at an 
arbitrary temperature. We have then derived the reduced density operator in closed form of the coupled oscillator in 
a non-equilibrium state at an arbitrary time. In doing so, we have applied both scalar-potential and vector-potential 
gauges for the interaction Hamiltonian. We believe that this exact expression for the reduced density operator will 
provide a useful starting point, e.g., for later useful discussions of quantum thermodynamics and quantum information 
theory within quantum Browian oscillator. 

The author would like to thank G. Mahler for helpful discussions. He also appreciates all constructive remarks of 
an anonymous referee. 
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Figure [l] y = {q(t))Y^ versus time t for a static field Fi{t) - Ai 0(f). The response {q(t))f^ is temperature- 
independent; refer to equations (43 i and (67 1. Here W() = Q = M=1. (1) dash (7 = 5, overdamped): From bottom 
to top, (black, Ai = 1), (violet, Ai = 1.5), and (red, Ai - 2); (2) solid (y = 0.1, underdamped): From bottom to top, 
(blue, Ai = 1), (maroon, Ai = 1.5), and (green, Ai = 2). 



Figure ^ 



y = {qit))f^ versus time f for an oscillatory field ^2(0 = A2 sinw/f. The response {q{t)}f^ is 
temperature-independent. Here wo = Q = M = A2 = 1; (1) violet dot: 7 - 0.1 (underdamped) and w/ = 1 ~ ojq, 
resonant [cf. ([65|)]; (2) blue dashdot: y = 0.1 (underdamped) and aij- - 1.5; (3) black dash: y = 5 (overdamped) and 
a>/ = 1; (4) red solid: y = 5 (overdamped) and cof = 1.5. 



Figure [3] y - {q^(t)}f^ versus time t for Fi{t) - Ai 0(r). For {q^it))f^ refer to equations (52ai and (63ai. 
Here h = ks - yvo - O. - M - 1. (1) black dashdot (y - 5, overdamped, and Ai = 1): From bottom to 
top, dimensionless temperature ksT/fvwo - 0.01,2,5; (2) blue dash (y - 0.1, underdamped, and Ai - 1): From 
bottom to top, ksT/fmo - 0.01,2,5; (3) red solid (y = 0.1, underdamped, and Ai = 2): From bottom to top, 
kBT/h-Wo = 0.01,2,5. 

Figure|4] y - {q^{t)}f^ versus time t for F2(t) - A2 sin w/f. Here h = kg = - - M = A2 = I. (l) black dot 
(y = 5, overdamped, and ojf - 1): From bottom to top, kBT/fiWQ - 0.01,2,5; (2) blue solid (y = 0.1, underdamped, 
and = 1, resonant): From bottom to top, ksT/frwo = 0.01,2,5; (3)reddash(y = 0.1, underdamped, and = 1.5): 
From bottom to top, ksT/hwo = 0.01, 2, 5. 



Figurejs] y - Dl(t) versus time t within the Drude damping model. For Z)j(f) refer to (68 1. Here h - kg - y^io — 
£1 = M = A] = A2 = 1, and y - 0.1, underdamped. (1) dash (for d.c. field): From top to bottom, dimensionless 
temperature ksT/fmo = 0.01,2, 5; (2) solid (for a.c. field): the same as for dash, with = 1 « wq, resonant. 
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Figure 2: 
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Figure 4: 
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Figure 5: 
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